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I. INTRODUCTION 

The problem of strong electron correlation is one of the 
most challenging problems in condensed matter physics. 
It is interesting because many materials with important 
properties for applications derive those properties from 
the delicate balance between minimizing the kinetic and 
the potential energies in strongly correlated materials. 
Most theoretical work on this problem has focused on 
equilibrium properties and linear response, with only lim- 
ited work available on the nonequilibrium system (which 
is most easily attained when driven by an external elec- 
tric field) . Two commonly studied models of strong elec- 
tron correlation are the Hubbard model 1 and the Falicov- 
Kimball models (equivalent to the Hubbard model with 
zero hopping parameter for the down spin electrons). De- 
spite tremendous efforts, the exact equilibrium solution 
of these models is known only in some limiting cases like 
one dimension, where the Bethe anzatz technique^ can 
be successfully applied to the Hubbard model, and in 
the infinite-dimensional case where both models can be 
solved^ with dynamical mean- field theory (DMFT). 

Since the exact solution of these problems is challeng- 
ing to attain, exact results in the form of sum rules can 
be quite valuable in determining the fidelity of different 
approximation techniques (be they analytic, variational, 
perturbative, or numerical approximations). This prob- 
lem was analyzed in equilibrium for the Hubbard model 
by Steve White& and used to check the accuracy of a 
quantum Monte Carlo solution to the two-dimensional 
Hubbard model. 

Can similar results be found for nonequilibrium sit- 
uations, like the case of a strongly correlated material 
in an uniform external electric field (with arbitrary time 
dependence)? The answer to this question has, surpris- 
ingly, not been discussed much in the literature. It is well 
known that the canonical anticommutation relation for 
fermion creation and annihilation operators leads to the 
integral of the spectral function being 1 in both equilib- 
rium and nonequilibrium situations. It is also well known 
that the conventional proof that the spectral function is 



nonnegative in equilibrium (arising from a Lehmann rep- 
resentation), does not apply to the nonequilibrium case, 
so the spectral function can be negative. It is also known 
that in the limit of a steady state, the spectral function 
recovers its nonnegativity. In this work, we examine the 
theoretical problem of the first few moments of the spec- 
tral functions in the presence of an external electric field. 
It turns out that the equilibrium results for the first few 
moments are quite similar to the nonequilibrium results, 
implying they can be used effectively to determine the 
accuracy of different approximation techniques in solv- 
ing nonequilibrium problems. 

Dynamical mean-field theory has been employed to 
solve many of the models of strongly correlated electrons. 
To date, most of this work has focused on equilibrium 
properties. Schmidt and Moment, made a first attempt 
to solve nonequilibrium DMFT via second-order pertur- 
bation theory for the Hubbard model. Their theoreti- 
cal formulation evaluated the case of a spatially uniform, 
but time dependent scalar potential, which unfortunately 
does not correspond to any electric field. More recently, 
a generalization of the Brandt-Urbanek solution^ for the 
localized electron spectral function allows for an exact 
numerical solution of the nonequilibrium problem for the 
Falicov-Kimball models within the DMFT framework. 
The procedure works directly in time by discretizing 
the continuous matrix operators and solving the nonlin- 
ear equations by iterating matrix operations on the dis- 
cretized operators. The approach has been tested against 
the equilibrium solution 9 !, and nonequilibrium results for 
the quenching of Bloch oscillations will be presented else- 
where. 

In this contribution, we derive operator identities for 
the first two spectral moments of the nonequilibrium 
Green's functions when the strongly correlated material 
is in the presence of a spatially uniform electric field with 
arbitrary time dependence. These identities are found 
for both the Falicov-Kimball and the Hubbard model. 
Our results are valid for any spatial dimensionality. In 
the general case, the moments depend on time, which 
is expected because the field can be turned on at any 
time; surprisingly, the first two moments of the local re- 
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tarded Green's function are time independent, and have 
the same form as in equilibrium. This last result is par- 
ticularly surprising because the local retarded Green's 
function is a nontrivial oscillating function of time in the 
noninteracting caseifi. 

The rest of the paper is organized as follows. We de- 
rive analytical expressions for the spectral moments for 
the Falicov-Kimball model and the Hubbard model in 
Section II. The application of nonequilibrium DMFT to 
the Falicov-Kimball model follows in Section III. Section 
IV contains results for the Green's functions and spectral 
moments in this case. A summary of our results and our 
conclusions appear in Section V. 



II. FORMALISM FOR THE SPECTRAL 
MOMENTS 



We begin our formal developments with a derivation 
of exact expressions for the zeroth, first and second spec- 
tral moments of the retarded and lesser Green's func- 
tions; the analysis is performed for the spinless Falicov- 
Kimball model and the spin one-half Hubbard model in 
an external arbitrary time-dependent spatially uniform 
electric field. The Hamiltonian for both models has the 
following form (in the absense of an external field) : 

no) = E'<^ i>y,r : <; 

ij ij i 



In the case of the Falicov-Kimball model, 2 the Hamil- 
tonian in Eq. describes a system which consists of 
two kinds of spinless electrons: itinerant c-electrons with 
a nearest neighbor hopping matrix ty and localized /- 

electrons with a hopping matrix equal to zero (t^ = 0). 
We normally take the hopping matrix to be between near- 
est neighbors only, but this is not a requirement. We 
do assume the matrix elements are all real and that the 
hopping matrices are Hermitian. The on-site interaction 
between the two electrons is equal to U. Double occu- 
pation by a c or / electron is forbidden by the Pauli 
exclusion principle. The chemical potentials are \x and 
/// for the c- and /-electrons, respectively. We will set 
(if = in our calculations; it plays no role in the spec- 
tral moments of the c particles, which we will be evalu- 
ating in this contribution. The spectral moments of the 
/-electrons in equilibrium were worked out in Ref. ITU 
In the case of the Hubbard modelji the Hamiltonian in 
Eq. describes a system of spin-up c = Cj-electrons, 
and spin-down / = Cj-electrons with equal hopping ma- 
trix elements t{j = t{j and chemical potentials /i = [if. 
The local Coulomb repulsion between electrons with dif- 
ferent spins is U. 

The electric field E(r, f) can be described by a vector 
potential A(r,t) in the Hamiltonian or temporal gauge 



(where the scalar potential vanishes): 



EM) 



1 dA(r, t) 

c dt 



(2) 



We assume that vector potential A(r, t) is smooth 
enough, that the magnetic field produced by A(r, t) can 
be neglected. For simplicity we shall assume that the 
vector potential is spatially uniform (independent of r). 
Note that if we have any time dependence to the elec- 
tric field, then neglecting the magnetic field will vio- 
late Maxwell's equations. In most situations these mag- 
netic field effects are small enough that they can be 
neglected in a first analysis, and added back later via 
either perturbative or gradient-based approaches. The 
electric field is coupled to the electrons via the Peierls 
substitution, 12 which involves modifying the hopping ma- 
trix elements by a phase that depends on the line integral 
of the vector potential: 
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The second line in each equation follows for spatially uni- 
form vector potentials. Note that the Hamiltonian in a 
field, Tt(A), is identical in form to that shown in Eq. QJ, 
but it uses the hopping matrices in Eqs. (|3J| and (J3J. Note 
also that v = for the Falicov-Kimball model. 

The "Peierls substituted" Hamiltonian in Eq. , with 
the hopping matrix elements in Eqs. Q and J2J, has a 
simple form in the momentum representation: 
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where the fermionic creation and annihilation operators 
now create or annihilate electrons with well-defined mo- 
mentum. The free electron energy spectra in Eq. @ 
satisfy 
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for the Hubbard model. In the case of the Falicov- 
Kimball model, the term vanishes. 

We shall consider the spectral moments for the re- 
tarded 

G£(*i,*a) = ~iO(h - ta) ({c*(h),cl(t 2 )}) (7) 
and the lesser 

G<(t 1 ,t 2 )=i(ct(t 2 )c k (t 1 )) (8) 

Green's functions; the symbol {0\, O2} = O1O2 + 2 0\ 
is the anticommutator and the operators cj^i) and Ck(i) 
are in the Heisenberg representation, where all time de- 
pendence is carried by the operators and the states are 
time-independent. Any Heisenberg representation oper- 
ator Or is connected with a corresponding Schrodinger 
representation operator Os via 



o H (t) = 



T exp<{ (i/h) I dtHi{i) 

to 



xc 



-a/h)H{o)(t-t ) 



e (i/K)H(G)(t-U>)Q s 
t 



T exp dtHi{t) 



,(9) 



where 7i(0) is the time-independent part of the Hamilto- 
nian [in Eq. with hopping matrix elements given by 
their field- free constant values], and Tti(t) is the time- 
dependent part of the Hamiltonian, which includes the 
interaction with an external field [but expressed in the 
interaction representation as detailed in Eq. I|1(J|I be- 
low]: that is, we define the time-dependent piece in 
the Schrodinger representation via His(t) = 'H(A) — 
Tt(0) and then re-express in the interaction representa- 
tion. Note that the interaction representation operator 
is defined to be the Schrodinger representation operator 
evolved under the time- independent Hamiltonian [i.e., 
the middle three terms in Eq. ©]. Hence, the time- 
dependent piece of the Hamiltonian in the interaction 
representation is expressed by 



Hi(t) 



o (i/h)H(0)(t-t ) 



Hi S (t)e 



-(i/h)H(0)(t-t o ) 



(10) 



in terms of the Schrodinger operator; in this form, there 
is the bare time dependence arising from the time depen- 
dence of the fields, plus the time dependence inherited by 
the operators, as we go from the Schrodinger represen- 
tation to the interaction representation. The symbol T 
(T) in Eq. @ is the time-ordering (anti-time-ordering) 
operator. 

We prepare our system to be in equilibrium prior to 
the field being turned on, hence the quantum statistical 
averages in Eqs. and ||SJ are defined with respect to 
the zero-field (equilibrium) Hamiltonian 7i(0): 



((...))= TV [ e -W°> (...)] /Z, 
where the partition function satisfies 
Z = Tr fe-^°) 



(11) 



(12) 



and P is the inverse temperature of the original equilib- 
rium distribution. As was already mentioned in Section 
I, the retarded and the lesser Green's functions form an 
independent Green's function basis. Any other Green's 
function can be expressed in terms of these two func- 
tions. This is in contrast to the equilibrium case, where 
only one Green's function is independent (because the 
Fermi-Dirac distribution function is determined by the 
equilibrium condition). 

Calculations in the nonequilibrium case are compli- 
cated by the fact that the Green's functions in Eqs. Q 
and JSJl are functions of two time variables, contrary to 
the equilibrium case, where they only depend on the rel- 
ative time difference (because the equilibrium system is 
time-translation invariant). It is convenient to transform 
the two-time dependence of the Green's functions from 
t\ and t 2 to Wigner coordinates, which use the average 
time T = (ti + t2)/2 and the relative time t = t% — t 2 
(do not confuse the average time T with the tempera- 
ture 1//3). Next, the relative time dependence is Fourier 
transformed to a frequency, and the additional (average) 
time evolution of different quantities is described by the 
average time coordinate T; in equilibrium, there is no T 
dependence. For example, the spectral function for the 
retarded and the lesser Green's functions can be defined 



as 



dte iuJt 



jjijlmG^Cr.f), (13) 
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where G^'' < (r, t) is the respective Green's function from 
Eq. or © with t\ = T + t/2 and t 2 = T - t/2; 77 is 
equal to —1 for the retarded Green's function and +1 for 
the lesser Green's function so that the spectral functions 
are nonnegative in equilibrium. In general, the spectral 
function depends on the average time, because the system 
no longer has time-translation invariance when a field is 
turned on at a specific time. 

We define the nth moment of the retarded and lesser 
spectral function [in Eq. (|13fl ] jJ.n y< to be 



dwwM^ ,< (T,w) 



J duio n fr?~) ImG^Cr,^) 



duju 11 77- 

7T 



/oo 
dte iwt Gk' < (T,t). (14) 
-00 

It is easy to show that this expression is equivalent to 
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Begin by noting that 

an 



so that 
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Substituting Eq. {HI into the first line of Eq. {H} then 
yields Eq. Q15(l. Before proceeding with the evaluation 
of analytical expressions for the spectral moments from 
Eq. 1151) . we note that the integration over frequency in 
Eq. HI 5(1 can be evaluated, yielding the following expres- 
sion, which connects the spectral moments to the deriva- 
tive of the Green's function with respect to relative time 
t at zero relative time: 



^• < (k,T) = r ? 2Im 



i n ^—G^' < (T,t) 



(18) 



This formula assumes that the Green's function is a dif- 
ferentiable function, which is true in most cases of inter- 
est. Despite the fact that the expressions in Eqs. 1 (15(1 
and 1(18(1 are formally equivalent, it is preferable to use 
one or the other in specific cases. 

In the case of the retarded Green's function, the well- 
known expression for the zeroth spectral moment can be 
found from Eqs. {TBJ) and Q: 



/i«(k,T) = -2Im [G£(T,t)] 
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= 1. 



T-\ 



(19) 



In the derivation of Eq. 1(19(1 . we used the anticommuta- 
tion relation for Heisenberg operators and the fact that 
the theta function is equal to 1/2 when its argument is 
equal to zero. 

It is more convenient to use Eqs. 1(15(1 and J7J to eval- 
uate the expression for the first moment of G R : 
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Taking the time derivative in Eq. 1(20(1 gives 
1 
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(21) 



The first term in Eq. ((21(1 is equal to zero, because the 
integral over time is equal to 1, therefore its imaginary 
part vanishes. The second term in Eq. 121(1 can be sim- 
plified by performing an integration over uj and replacing 
the time derivatives of the operators by their commuta- 
tors with the Hamiltonian, according to the Heisenberg 
equation of motion idO(t)/dt = [0(t), 'Hit)], where 7i(t) 
is the total Hamiltonian including the effects of the time- 
dependent field. This yields 

M ?(k,T) = iRe(({[ck(T),W(T)],4(T)}) 

- ({c k (T),[4(T),W(T)]})). (22) 

Evaluation of the commutators of the Fermi-operators 
with the Hamiltonian and the subsequent anticommuta- 
tors in Eq. J^J gives the following expression for the first 
spectral moment of the retarded Green's function 



//f(k,T)=e k 



eA(T) 
he 



fx + Urif, 



where 



/ = £</k(nfk(T)) 



k 



(23) 



(24) 



is the average number of / (cjj-electrons in the system; 
this number of electrons does not depend on the average 
or the relative time, because the total electron number 
for each species of electron is conserved. 

Similarly, the expression for the second moment of the 
retarded Green's function can be found from Eqs. 1(15(1 
and J7J): 

»*(k,T) = iRc(({[[ Ck (T),H(T)],H(r)],ct(r)}) 

- 2({[c k (T),W(T)],[4(T),W(T)]}) 

+ ({c k (T),[[4(T),W(T)],W(T)]}}). (25) 

Details of the derivation are presented in the Appendix. 

Evaluating the commutators and anticommutators in 
Eq. gives 



M?(k,T) = 



+ 2U 



eA(T) 

he 
he 
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(26) 



nj + U 2 rif. 



The moments of the local retarded Green's function 
fi R (T) are obtained by summing the corresponding spec- 
tral moment functions /U*(k, T) over k 

£*(T)=5>*(k,T). (27) 

k 

Performing the summations for Eqs. 1(23(1 and ((26(1 yields 
the following local moments: 

A?(T) = -fi+Unt; (28) 
1 



+ /r - 2U(xn,f + U z n f . 



(29) 
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These results coincide with those derived previously for 
the Hubbard model in the equilibrium c&sem. Since the 
hopping matrix is always chosen to be traceless in our 
models, the sum of the energy ^ k e(k) in Eqs. ill'- ill and 
(|26|l is equal to zero. The expression for the zeroth lo- 
cal moment fiy has the same form as the expression for 
the zeroth spectral moment in Eq. I|19fl . since the zeroth 
spectral moment is momentum-independent. Hence, the 
zeroth and the first two local moments of the retarded 
Green's function in an arbitrary external time-dependent 
homogeneous electric field are all time-independent! This 
is a nontrivial result because the retarded Green's func- 
tion strongly depends on the average time. In particular, 
the retarded Green's function is an oscillating function 
of timeifi when [7 = 0. Furthermore, the moments do 
not depend on the electric field at half-filling because the 
chemical potential is not changed by the field! It isn't ob- 
vious whether the chemical potential would be changed 
by the field off of half-filling. 

In the case of half- filling, where n c — rif = 1/2 and 
^ = £7/2, the expressions in Eqs. JTHJ), J2E1) and (f2l3|) 
acquire an even simpler form: 



£f(T) = 



1, 
0, 

1 u 2 

2 + - 



(30) 
(31) 

(32) 



If one examines the moments for gauge-invariant 
Green's functions^ then the local moments are un- 
changed, and the spectral function moments are modified 
by a time-dependent shift of the momentum wavevector. 
We don't include those formulas here, because they just 
involve such a simple shift. 



The corresponding moments of the lesser Green's func- 
tions are found by a similar analysis. Using Eqs. (|18f) and 
© we find 



/x<(k,T) = 2n c (k,T), 
M <(k,T) = -Re(([4(T),W(T)] Ck (T) 

[ci(T)[c k (T),n(T)} 



(33) 



(34) 

l4(k,T) = iRe(([[ct(T),W(T)],H(T)] Ck (T)) 
- 2([4(T),W(T)][c k (T),W(T)]) 
+ U(T)[{ct(T),H(T)},H(T)})) , (35) 



where 



n c (k,T) = (4(T)c k (T) 



(36) 



is the momentum distribution function for the c (c-|-)- 
electrons. Note that a commutator term depending on 
the derivative of the Hamiltonian with respect to time 
can be shown to cancel, so it is not included in the second 
moment expression above. Evaluation of the commuta- 
tors of the operators with the Hamiltonian [in Eq. (JjJJ] in 
Eqs. (EH) and gives 
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In Eq. I|38|) . we have suppressed the time label T corre- sponding to the time at which all operators are evaluated. 
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We continue to suppress this time label in some equations 
below; this should not cause any confusion. 

The expressions in Eqs. (|33[l - l|38[l for the lesser spectral 
moments are more complicated than the corresponding 
retarded moments. However, they simplify in the case of 
the local Green's function, where we find 
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k 



(k,T) 



(41) 



(42) 



being the time-independent particle density of the c (c-|-)- 
electrons. In order to save space, and make the equations 
more transparent, we use a mixed real-space/momentum- 
space representation for the operators in Eqs. I|40|) and 
l)4ip. Note that the first moment involves one correlation 
function and the second moment involves two correla- 
tion functions. Note further that the value of the first 
moment in Eq. (|40l) is equal to twice the average value 
of the Hamiltonian. The last term in Eq. (|41|) is equal 
to zero in the case of half- filling (/i = £7/2). The second 
two moments of the lesser Green's functions do appear to 
depend both on the average time and the electric field (al- 
though our empirical evidence at half filling suggests the 
second moment may be independent of average time — see 
the numerical results below). 

There is an interesting observation that can be made 
about the first moment, and its relation to the current 
driven by the electric field and the phenomenon of Bloch 
oscillations. In the limit where U is small, one can eval- 
uate the correlation function in Eq. (|40|l via a mean-field 
theory decoupling ((///jctc,) w (///,) (cjcj)). For ex- 
ample, at half filling, the first moment will be equal to 
twice the average kinetic energy (including the shift by 
the vector potential needed to construct the actual ki- 
netic energy from the bandstructure) plus a correction of 
order U 2 . If the current oscillates, we expect the aver- 
age kinetic energy to oscillate as well. Hence, for small 
U there is a correlation between oscillations in the first 
moment of the local Green's function and oscillations of 
the current. 



The correlation function i^f} fic\c-ij that appears in 

Eqs. (|40|l and (|41|l can be determined for the Falicov- 
Kimball or Hubbard model via the equation of motion, 
because it is related to the total energy of the Hamil- 
tonian. To show how this works, we first provide the 
derivation for the equilibrium case using an imaginary- 
time formalism. Begin with the definition of the Green's 
function in real space 



Gijir) = - (T Ci (T)c](0)) , 



(43) 



with a similar result for the spin-down electrons 
in the Hubbard model. Here we have Cj(r) = 
exp(7i-r)ci(0) exp(— Tir). Taking the imaginary-time 
derivative of the local Green's function gives 

d T Gu(r) = -S(r)-(T[H,Ci(T)}c\(0)) 

= -8{r) + t^Gi + si{T) + nGii{T) 
s 

+ u(TfJfMr)cl(0)), (44) 



where the symbol 6 denotes the translation vector to a 
nearest- neighbor site and i + S as a subscript refers to the 
lattice site that is the nearest neighbor of site i translated 
by the nearest-neighbor translation vector S. Hence we 
determine the correlation function via 



U(Tf}f4a) = lim [-d T G u (T) 



(45) 



tY,G i+Si {r) + fiG u (T) - 5(r) 



Using the Matsubara frequency representation 



G(iuj v 



Gn 







dre 1 



r G{r), 



(46) 



with iu> n — m(2n + l)//3 the fermionic Matsubara 
frequency, allows us to determine a simple expres- 
sion for the correlation function. Note that G(r) = 
(1//3) J2 n exp(— iuj n T)G n and that the Green's function 
in momentum space satisfies G„(k) = \/[iuj n + fi — 
S„(k) — e(k)], to find that the correlation function sim- 
plifies to 



S„(k) 



/i-E n (k)-e(k)' 



(47) 

In the limit of infinite dimensions, the self-energy is a lo- 
cal function, and hence has no momentum dependence. 
Then the sum over momentum can be performed by 
changing from a sum over momentum to an integral over 
the noninteracting density of states. This produces the 
local Green's function, and we are left with the final form 
for DMFT: 



1 



(48) 
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In numerical calculations, it is more convenient to eval- 
uate the summation in Eq. 148f) via the formally equiva- 
lent expression with the Hartree-Fock contribution to the 
self-energy removed 



Tftfic\ Ci 



:/V>v-4 + ^£[^-tf</V>]G n , 



(49) 

because the Matsubara summation converges faster. 

In the noncquilibrium case, one can perform a similar 
analysis, but now one has to work with the nonequilib- 
rium Green's functions as functions of real time variables. 
Using the standard equations of motion, and definitions 
for nonequilibrium Green's functions, one can show, af- 
ter some significant algebra, that the correlation function 
can be expressed as 



,d_ 
''dtl 



eAjtx) 
he 



x G<(h,t 2 ) 



iJ2 f dt[l$(h,t)G2(t,ti) 

1, J 



+ £f(ii,t)G£(t,ti)], 



(50) 



where the final result is written in terms of retarded, 
lesser and advanced Green's functions and self-energies 
in the presence of the field. In the DMFT limit, the self- 
energies have no momentum dependence. This expres- 
sion appears like it can have average time dependence, 
but we cannot say that it definitely does, because there 
could be a cancellation of the time dependence. 

As a check of Eq. (|50|l , we evaluate it in equilibrium, to 
show it yields the same result as Eq. (|4^|) . When we are 
in equilibrium, the correlation function is independent of 
time and the Green's functions and self-energies depend 
only on the time difference of their two time variables. 
Hence, we can perform a Fourier transform by using the 
convolution theorem to transform Eq. (|50|l into 



U(fjfi4c, 



i 

2^ 



J2 J <m s £m^m 



(51) 



Using the fact that the lesser functions satisfy G^ = 
-2if(u>)1mGg(u) and £< = -2if(w)hnEg(u) in equi- 
librium, allows us to transform Eq. (|51|l into 

U(f!fi4ci) = ~J2 I duf(u)Im [Ef(w)G£M] , 

^ k J 

(52) 

which is equal to the analytic continuation of Eq. I|48|) 
from the imaginary axis to the real axis. 

The correlation function ^/p+ q /pCjj:_qCky which en- 
ters Eq. 141|) is more complicated to evaluate, because 



it cannot be expressed in terms of a simple equation 
of motion for the single-particle Green's functions. Be- 
cause of it's complex nature, we will evaluate it only for 
the Falicov-Kimball model in equilibrium. While it is 
certainly true that it can be evaluated for the Falicov- 
Kimball model in the presence of a field, some of the for- 
mal details become quite complicated, and take us away 
from the main theme of this work, so we do not perform 
such an analysis here. 

Our starting point, then, is the operator average 



£ [e(k) + e(k-q)](/ P+q / pC t_ q c k 



k,p,q 



(53) 



where we have set the vector potential A equal to zero. 
We use a Fourier transform to express the localized 
electrons in terms of their real-space operators. Then 
Elitzur's theorem^ ensures that the operator expectation 
value vanishes if the two localized electrons are not at the 
same lattice site (i.e., there is no spontaneous hybridiza- 
tion in the Falicov-Kimball model for nonzero tempera- 
ture). This allows us to perform the summation over the 
momentum variable p and gives us 



]T ]T [e(k) + e(k - q)] (///4_ q c k 

k,q i 



(54) 



Next, we express the bandstructure in terms of the 
summation over nearest neighbor translation vectors 5: 
e(k) = — t* ^2 S exp[ik ■ 6]/Vd, and introduce Fourier 
transforms for the itinerant electrons to real space. This 
allows us to sum over the remaining momenta, yielding 



f f f c f 



i+5 1 



(55) 



The statistical averages in Eq. (|55|l have already been 
evaluated^ The procedure is to imagine adding a small 
field — Yli hifjfi to the Hamiltonian, and evaluate the 
expectation value with the localized particle number via 
a derivative with respect to the field strength hi (then 
set the field to zero to evaluate the average). This gives 



—y 



iS 



— T 

id 

[G l+ si(T 



1 d 

Pdhi 

1 d 

-> 0-) + G 11+s (t 



CiCi +S ) + (C [ i+s Ci 



{Wi 



Now we follow the derivation in Ref. [l 
the derivatives from the following: 



l_d_ 

]3dh. 



G 1+Si {t^Q-) 



n jk 

-y 



1 d 



1 d 



Gi+si -r — S(iw n )Gji (iui ri 
p oh 



{)-)] . (56) 
^, which evaluates 

(57) 

G~^{iuj n )G kl {iuj n ) 



U ] f)G u {iu n ) 
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where the derivative with respect to the field h acts on 
the local self-energy. After some long and complicated 
algebra, that derivative can be determined, which yields 
our final result for the correlation function 



I mi ^Gk(w n )£(wj„) 



U 



n k 



pu 



5^[-l + (tw n + M-E n )G n ]E n . (58) 



This result is quite similar to the previous result for the 
other correlation function, except now we have an extra 
weighting factor of 2e(k) in the summation over momen- 
tum. 



III. THE FALICOV-KIMBALL MODEL IN 
INFINITE DIMENSIONS 

In this Section, we examine the time-dependence of the 
local Green's function for the Falicov-Kimball model on 
an infinite-dimensional hypercubic lattice. We consider 
the case of half-filling with the system being coupled to 
an external homogeneous time-dependent electric field. 
The time dependence is taken to be particularly sim- 
ple, at t — to a constant field is instantly turned on. 
The formalism involves generalizing the DMFT to the 
nonequilibrium case. The way to do this is based on 
a Kadanoff-Baym approach in real time, where continu- 
ous matrix operators are discretized along the Kadanoff- 
Baym contour, and operator manipulations are carried 
out on the discretized matrices using standard linear- 
algebra approaches. A short description of this tech- 
nique, including a benchmark against the well-known 
equilibrium solutions has already appeared^; further de- 
tails of this approach will appear elsewhere. To test our 
formulas for the first few moments, we compare numeri- 
cally calculated local moments to the exact moments in 
Eqs. (|5UP21) and ElgTJ- 

The action for the Falicov-Kimball model is quadratic 
in the conduction electrons. Hence the Feynman path in- 
tegral over the Kadanoff-Baym contour can be expressed 
by a determinant of a continuous matrix operator whose 
arguments are defined on the contour. Since the con- 
centration of static particles on each site is conserved, 
the trace over the fermionic variables can be straight- 
forwardly taken. This is what allows the nonequilib- 
rium DMFT problem to be solved, but the technical de- 
tails are complicated. The first thing that needs to be 
noted is that the self-energy remains local even in the 
presence of a field. This follows by applying Langrcth's 
rulesi^ to the perturbation theory, which state that every 
nonequilibrium diagram can be related to an equilibrium 
diagram, but now one must perform the analysis over 
the Kadanoff-Baym contour rather than over the finite 
imaginary time interval. Since the perturbative analysis 
of the equilibrium problem shows the self-energy to be 



loca in equilibrium, it remains local in the nonequi- 

librium case as well. 

We couple the system in Eq. JSJ to an external electric 
field along the unit-cell diagonal direction in real space; 
this yields the following vector potential for the electric 
ficldia: 



A(t) = A(t)(l, !,...,!). 



(59) 



The bandstructure for noninteracting electrons coupled 
to the electric field in Eq. (|59|) has a simple form: 



eA(t) 

he 



he J 



he J 



with the energy functions defined to be 
t* 



and 



£k 



-—= cos(ako) 



,.19 



tk- 

(60) 
(61) 

(62) 
i the 



and t* is the renormalized hopping parametei" 
limit of d — ► oo; we take t = t*/2Vd. 

Because the Green's functions now depend on two en- 
ergies, the summation over the infinite-dimensional Bril- 
louin zone can be replaced by a double integral over a 
joint density of states for the two energies. This joint 
density of states (DOS) isi 



P2{e,e) 



1 



irt* 2 a d 



exp 



£ 



£ 

^2 



(63) 



The numerical integration over the joint DOS is per- 
formed by an averaged Gaussian integration with 54 and 
55 points for each energy axis: 



cfeexp(— e 2 )F{e) 



N 

E 



w i F{e i ) 1 



(64) 



where Wi are Gaussian weights which correspond to the 
N energy points £j. Since the Green's functions often 
depend on the energy as exp(?c£), the Gaussian quadra- 
ture rule in Eq. H64fl fails to give correct results when c 
is on the order of (or larger than) the inverse of the grid 
spacing of the energy points near e = 0. In this case, 
the sum over discrete points contains a systematic con- 
tribution of terms which do not cancel each other, and 
leads to an overestimated value to the integral. One way 
to efficiently correct this is to average two Gaussian sum- 
mations with numbers of Gaussian points equal to N and 
iV + 1 because the Gaussian points interleave each other, 
and act like a step size about half as big as either sum 
alone, and they give somewhat better accuracy than per- 
forming the integral with 2N+ 1 points, because a subset 
of those points are at such large absolute value that the 



9 



Gaussian weight is small enough that it can be safely 
neglected. 

In order to calculate the nonequilibrium local Green's 
functions, one needs to self-consistcntly solve a system of 
equations^ which connects these functions with the cor- 
responding local self-energy £(ti,t2) and an effective dy- 
namical mean-field X(ti,t2)', these equations are similar 
in form to the equilibrium case^*^ but now all the func- 
tions are discrete matrices of two time variables defined 
on the Kadanoff-Baym contour (see Fig. Details of 
the algorithm and the nonequilibrium DMFT equations 
will be given elsewhere^. 



'max 







'max 














"'max 'P 



FIG. 1: Kadanoff-Baym contour for the two-time Green's 
functions in the nonequilibrium case. We take the contour to 
run from —t max to t max and back, and then extends down- 
ward parallel to the imaginary axis for a distance of j3. The 
field is usually turned on at t = 0; i.e., the vector potential is 
nonzero only for positive times. 

First, we present results for the local moments in equi- 
librium, when the system of the DMFT equations is 
solved in the frequency representation using the Brandt- 
Mielsch approach^ (for details, see Ref. 0) . Plots of the 
local retarded and lesser DOS for different values of U are 
shown in Figs. [3 and [3] The metal-insulator transition 
occurs at U = \/2; the insulator has anomalous prop- 
erties because there is no real gap — instead the DOS is 
exponentially small in a gap region around the chemical 
potential, and vanishes only at the chemical potential. 

The moments for the retarded and lesser Green's func- 
tions are calculated by directly integrating the Green's 
functions multiplied by the corresponding power of fre- 
quency; we use a step size of Aw = 0.001 and a rect- 
angular quadrature rule. The results for half-filling are 
presented in Table [I] One can immediately see that 
the numerical results for the first moment of the lesser 
Green's function are in an excellent agreement with the 
exact expressions for the moment (the zeroth and sec- 
ond moments agreed exactly with their exact expres- 
sions). We also calculated the retarded moments and 
they agreed exactly with the exact values, so we don't 
summarize them in a table. Note that there appears 
to be a relation between the second moments of the re- 
tarded and lesser Green's functions at half filling. This 
relation ceases to hold off of half filling. For example, 
in the case with U = 2, w% = 0.25, and p e = 0.75 we 
find the following moments (exact results in parentheses): 
/if = -0.872072144 (-0.872071963); ^ = 2.01050948 




Frequency co [t*] 

FIG. 2: DOS for the equilibrium retarded Green's function for 
different values of U. The DOS is independent of temperature. 
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FIG. 3: Local lesser Green's function in equilibrium at /3 = 10 
and for different values of U. 



(2.01050951); fif = -2.149496 (-2.149462); and jlf = 
3.623019 (3.622970). These results are more indicative of 
the general case (where ^ jlf). Note that one needs 
to use many Matsubara frequencies in the summations 
to get good convergence for the average kinetic energy 
and for the second correlation function (we used 50,000 
in this calculation with /3 = 10). The majority of our 
numerical error comes from the difficulty in exactly cal- 
culating those results; indeed, the exact result, calculated 
from the operator averages on the Matsubara frequency 
axis is probably less accurate than the direct integration 
of the moment on the real axis. One can improve this 
situation somewhat by working on the real axis to cal- 
culate the different operator averages, but we wanted to 
indicate the accuracy under the most challenging circum- 
stances. We feel our final results are quite satisfactory, 
and indicate these sum rules do hold. 

We now focus on numerical results for the nonequi- 
librium code. Our first benchmark is to calculate equi- 
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TABLE I: First spectral moment for the lesser Green's function in equilibrium with /3 = 10 and different values of U. The 
zeroth moment is accurate to more than eight digits, and is not included. Similarly, we find the second moment is equal to 
0.5 + U 2 /4 to high accuracy, and is not included. Note that the first moment continuously evolves from the value — 1/i/tt for 
U = and j3 — oo to approximately —(7/2 as U increases. The approach to —U/2 is expected due to the formation of upper 
and lower Hubbard bands separated by U. 



moment 


U = 0.5 


U = 1.0 


U = 1.5 


U = 2.0 


U = 1.0 


U = 6.0 


P? 


-0.591699 


-0.717901 


-0.902869 


-1.119047 
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-3.041526 


fif (exact) 


-0.591687 


-0.717886 
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FIG. 4: (Panel a) Imaginary part of the lesser Green's func- 
tion as a function of the relative time coordinate for different 
discretizations of the time contour. The model parameters 
are U = 0.5, f3 = 10 and E — 0. The average time T is set 
equal to zero. The parameters for the Kadanoff-Baym time- 
contour discretization are: tmax = 15 and At = 0.1 (i.e., 
100 points taken along the imaginary axis); the discretization 
along the real time axis is given by At as shown in the figure. 
Note how the results systematically approach the exact result 
as the discretization goes to zero. (Panel b) Imaginary part 
of the lesser Green's function as a function of frequency for 
different discretizations of the time contour. 



librium results with that code and compare with ex- 
act results available for the equilibrium case. Such 
calculations^ show good convergence and precision when 
U lies below the critical U for the metal-insulator tran- 
sition (U < v2)' Here we demonstrate this by showing 
how the relative time dependence of the imaginary part 



TABLE II: Spectral moments for the retarded Green's func- 
tion in the case of zero electric field at U — 0.5, /3 — 10 and 
calculated with the nonequilibrium real-time formalism with 
different values for the time step At. The other parameters 
are t max =15, N = 54, 55, At = 0.1 



moment At = 0.1 At = 0.067 At = 0.05 At = 0.033 exact 
1.580785 1.331640 1.232022 1.144811 1 
0.174040 0.082610 0.052785 0.030002 
1.324976 0.979230 0.848047 0.737020 0.5625 



Ao 
A? 

PS 



TABLE III: The same as in Table III but for the case of the 
lesser Green's function. 



moment 
Ao 
Pi 
P2 



At = 0.1 At = 0.067 At = 0.05 At = 0.033 exact 
1.480893 1.289036 1.207662 1.133850 1 
-1.036753 -0.850675 -0.774525 -0.706893 -0.591687 
1.108705 0.870853 0.777152 0.695791 0.5625 



of the lesser Green's function at U = 0.5 and j3 = 10 
changes when one decreases the time step At on the real 
part of the Kadanoff-Baym contour. As follows from 
Fig. 01 the solution becomes more accurate as At de- 
creases but the accuracy is reduced as the temperature is 
lowered, as can be seen by comparing to the = 1 results 
in Ref.|9l For this calculation, we determine the moments 
by taking the derivatives of the Green's functions in the 
time representation. The values for the spectral moments 
at [7 = 0. 5, /3 = 10 and different values of At are pre- 
sented in Tables ITT1 and Hill The results for the moments 
improve as At decreases. 

In the insulating phase, the calculations are less accu- 
rate (see Fig. |SJ| . The self-energy develops a pole in the 
frequency representation, which gives the imaginary part 
a delta function at the chemical potential. Hence, we ex- 
pect a constant background value for the self-energy in 
the time representation. This cannot be properly repre- 
sented in a numerical calculation that has a finite cutoff 
in the time domain, which makes the real-time formal- 
ism much more challenging in the insulating phase. The 
problems appear to be somewhat less pronounced for the 
Green's functions, but the convergence is slower than in 
the metallic phase, and the gap region has unphysical 
behavior [the DOS goes negative in the gap region as 
shown in the inset to panel (b)]. One cannot get rid of 
these oscillations without having a time-domain cutoff 
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FIG. 5: (Panel a) Imaginary part of the lesser Green's func- 
tion as a function of the relative time coordinate for different 
discretizations of the time contour. The model parameters 
are U = 2, j3 = 10 and E = 0. The average time T is set 
equal to zero. The parameters for the Kadanoff-Baym time- 
contour discretization are: tmax = 15 and At = 0.1 (i.e., 
100 points taken along the imaginary axis); the discretization 
along the real time axis is given by At as shown in the figure. 
Note how the results systematically approach the exact result 
as the discretization goes to zero. (Panel b) Imaginary part 
of the lesser Green's function as a function of frequency for 
different discretizations of the time contour; in the inset, the 
region around u = is blown up to show the gap develop- 
ment as a function of the discretization. Note how the gap 
region converges very slowly — instead we see the DOS go neg- 
ative in the gap region. Properly determining that structure 
in the frequency domain requires the Green's function over 
an extended time domain, which is not numerically feasible. 



that extends to infinity, but by comparing with the ex- 
act results, we find that a time-domain cutoff of about 
t re i w 200 provides quite reasonable results for the calcu- 
lations; in our results with the nonequilibrium code, we 
are limited to time-domain cutoffs of closer to 15 — 30, 
which explains the poor agreement for the gap region. 
Fortunately, these numerical problems appear to reduce 
when an external field is turned on, and we are in the 
nonequilibrium case. 

Note that we show equilibrium results only for the av- 
erage time T = 0. In equilibrium, the results should be 



independent of T, but we find that we have a modest 
T dependence due to discretization error. The results at 
T = turn out to be the least accurate, and the accuracy 
of the results improves as the discretization size is made 
smaller. In general, we find the T dependence of the 
Green's functions to vary (pointwise) by no more than 
40% for At — 0.1 and to be reduced to a 10% variation 
when At = 0.033 (for U = 0.5). The variation is about 
three times larger for U = 2. We find less variation in 
the nonequilibrium calculations, which appear to be bet- 
ter suited to the real-time formalism [most likely because 
the Green's functions don't behave like complex expo- 
nentials exp(iei) as the equilibrium functions do; such 
functions can be particularly difficult to deal with in our 
real-time numerical calculations] . 

As a nonequilibrium problem, we consider the case 
of the interacting Falicov-Kimball model in the metallic 
phase at U = and U = 0.5 with /3 = 10 and a constant 
electric field E — 1 is turned on at some moment of time. 
The Green's functions for noninteracting electrons^ [see 
panel (a) of Fig. [S] are oscillatory functions of the rel- 
ative time coordinate in presence of the external field. 
The average time dependence is weak for relative times 
up to about 2T, after which, the Green's function decays 
as a function of t re \. Notably, the results for different 
average times lie on top of each other until the relative 
time becomes larger than 2T . This implies that there 
is little or no average time dependence to the retarded 
Green's function as T — > oo, and the retarded Green's 
function becomes a periodic function of t re i (with the 
Bloch period); this latter result implies that the Fourier 
transform will consist of evenly spaced delta functions, 
which is the familiar Wannier-Stark ladder. When inter- 
actions are turned on, the situation changes [see panel 
(b) of Fig.E]. We still see little average time dependence 
for relative times smaller than 2T, but the Green's func- 
tion does not appear periodic in t re i for small times. The 
critical question to ask is, what happens as T — > oo? If 
the Green's function becomes periodic in t re i , then there 
will be delta functions surviving in the Fourier transform, 
but if it continues to decay, there will not. Our data does 
not extend far enough out in time for us to be able to 
resolve this issue. 

It is also interesting to examine the density of states 
in frequency space. The U = case has been studied 
exhaustively^, so we don't repeat it here. The interact- 
ing case is plotted in Fig. [7] In constructing this plot, 
we set the real part of the retarded Green's function 
to zero before performing the Fourier transform, since 
the exact result must vanish by particle-hole symmetry, 
and our numerical calculations have a small nonzero real 
part. Note how the DOS rapidly readjusts itself into a 
nonequilibrium form, and how the steady state appears 
to be approximately reached. The only subtle issue is the 
one discussed above, of whether there will be delta func- 
tion peaks emerging in the DOS as the average time gets 
larger. What is clear, is that if such peaks do form, they 
are quite unlikely to appear at multiples of the Bloch 
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FIG. 6: Imaginary part of the retarded Green's function as a 
function of the relative time coordinate at different values of 
the average time (starting at T — and running to T = 40 
in steps of 5) . The field is switched on at the time T = and 
we take At = 0.1 along the Kadanoff-Baym contour. The 
model parameters are f3 = 10 and E = 1. Panel (a) shows 
the noninteracting result U — 0, while panel (b) shows the 
interacting result U = 0.5 (note that the curves extend as 
far out in t re i as we have data; the cutoff in t re i comes from 
the finite time domains of our calculations). Note how the 
results appear to retrace themselves for different average times 
until the relative time becomes larger than approximately 2T, 
where the Green's function decays. For the noninteracting 
case, the pattern obviously becomes periodic in the Bloch 
period as T — > oo, which leads to delta functions in the Fourier 
transform (with respect to t re i), whereas the interacting case 
may or may not be approaching a periodic form for large 
relative time; the data does not extend far enough out to be 
able to make a conclusion about the asymptotic form. 



frequencies, because we see no sharp peaks forming near 
integer frequencies here (the Bloch frequency is equal to 
1 for E = 1); indeed, the DOS seems to be suppressed at 
integer frequencies. 

The lesser Green's functions are remarkably similar to 
the retarded Green's functions, except they are nonzero 
for all t re i, not just for positive values. There also is 
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FIG. 7: Density of states for different average times (run- 
ning from T — to T = 40 in steps of 5) for the interacting 
system, with U = 0.5 and E = 1. These results correspond 
to the Fourier transform of the data in Fig. [5] Note how 
the DOS seems to be approaching a limiting form even for 
this small a value of the average time. We cannot tell, how- 
ever, whether there might be some low-weight delta functions 
appearing somewhere in the spectrum. In any case, it does 
appear that there are no sharp structures near the Bloch fre- 
quencies, which are at integer frequencies for E = 1. 



limited average time dependence, except in the region of 
small t re i , where the first derivative of the Green's func- 
tion does vary with average time. It is this variation 
that leads to oscillations in the current, and is critical for 
understanding the behavior of these systems. Unfortu- 
nately, it is not simple to directly evaluate such deriva- 
tives accurately when we calculate them numerically, be- 
cause our time step is rather large. We examine them in 
this numerical fashion, with results plotted in Fig. |H1 for 
the case At = 0.1. Note how the first moment oscillates, 
then decays, and finally seems to reach a steady oscil- 
latory state as the average time increases. We cannot 
tell whether the moment becomes a constant at large av- 
erage times or continues to oscillate from the data that 
we have, although it appears to be approaching an os- 
cillatory steady state. Note further that because U is 
not too large here, we anticipate a correlation between 
these oscillations in the first moment and oscillations of 
the current. Indeed, if one calculates the current, one 
finds that it also appears to approach a steady oscillat- 
ing state for large average time; details of these results 
will be presented elsewhere. 

Despite a significant change in the Green's functions 
after switching on an electric field (compare the T = 
results, which are close to the equilibrium results to the 
larger T values in Figs. El and UJ, the spectral moments 
for these functions do not change much (Tables IIVI and 
rVjl. The spectral moments are connected with the rela- 
tive time derivatives of the Green's functions [Eq. jlSjl]: 
the fact that some moments do not change in the pres- 
ence of an electric field suggests that the those derivatives 
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FIG. 8: First moment of the local lesser function jlf plotted 
as a function of average time T for U = 0.5 and E = 1. The 
field is turned on at T = 0. The step size is At = 0.1, and 
the moment is calculated from the numerical derivative of the 
data. 



TABLE IV: Spectral moments for the retarded Green's func- 
tion in the case when the constant external electric field E — 1 
is switched on at T = 0. The parameters are U — 0.5, /3 = 10, 
At = 0.05, tmax = 15, N = 54, 55, and At = 0.05. These 
moments should all be independent of time, and they appear 
to be within the numerical error. 



moment 


T=0 


T=5 


T=10 


T=15 


T=20 


exact 


fro 


1.0025 


1.0088 


0.9985 


0.9951 


0.9967 


1 


~R 

I'i 


0.00665 


-0.00054 


0.00005 


0.00054 


0.00003 





IJ-2 


0.56155 


0.56198 


0.55184 


0.55030 


0.55112 


0.5625 



are independent of the electric field. Indeed, we calcu- 
late the moments in these examples from the derivatives 
of the Green's functions at t = because we often do 
not have data out to a large enough relative time to per- 
form the Fourier transform, and evaluate the moment in 
the conventional way. Note that, the second moment, 
which is actually equal to the curvature of G R and G < , 
is independent of average time for the retarded function, 
and does not appear to have average time dependence for 
the lesser function. The first moment of the lesser func- 
tion does depend on average time (see Fig. |5J). Finally, 
it is interesting to observe that the values for the mo- 
ments are much closer to the exact results for a similar 
discretization size than what we found in the equilibrium 
case. This is why we believe that the real-time numerical 
algorithm converges better for the nonequilibrium case 
than for the equilibrium case. 

To conclude our numerical analysis, we consider the 
spectral moments of the lesser Green's function in 
the case where it is approximated by the generalized 
Kadanoff-Baym (GKB) ansataSi. The idea of the GKB 
is to represent the lesser Green's function in terms of 
the distribution function (determined by the t — limit 
of the lesser Green's function) and the full two-time re- 



TABLE V: The same as in Table IV for the case of the lesser 
Green's function. The first moment appears to change with 
average time in the nonequilibrium case. The equilibrium 
value is —0.5917, which agrees well with our result before 
the field is turned on. The second moment appears to be 
independent of average time, even in a held, but we have no 
proof of this fact. 



moment 


T=0 


T=5 


T=10 


T=15 


T=20 


exact 


Ao 


1.0025 


1.0098 


0.9997 


0.9960 


0.9975 


1 


t<l 


-0.5878 


0.0042 


-0.1618 


-0.0454 


0.0933 


? 


i4 


0.5520 


0.5666 


0.5554 


0.5572 


0.5588 


? 



tarded and advanced Green's functions: 

-i [G R (t u t 2 )G<(t 2: t 2 ) 
G<{h,t x )Gt{tx,h)} ■ (65) 



G£(fi,* 2 ) 



The expression for the spectral moments of the lesser 
Green's function, approximated by Eq. (|65|) . can be found 
by taking relative time derivatives of G k (T, t) as shown in 
Eq. i|18|l . Similar to the derivation performed in Section 
II, one finds (the hat denotes the GKB approximation) 



Ao ( k > T ) 
/}<(k,T) 



2n c (k,T), (66) 
n c (k,T) (({[ Ck (T),tf],4(T)}) 

{c k (T),[4(T),tf]} 
dn c (k,T) 



Im 



dT 



(67) 



A 2 < (k,T) 



±n c (k,T) (({[[c k (T),H],H],cl(T)} 

2({[c k (T),i/],[c k (T), J ff]} 

{c k (T),[[cl(T),H},H}}) 
d 2 n c (k,T) 



-Re- 



ar 2 



(68) 



Comparison of the expressions in Eqs. (|66[) -H68 [I for the 
GKB moments with the corresponding exact expressions 
in Eqs. {TH1), an d I23> f° r the moments of G R allows 
us to connect the lesser GKB spectral moments with the 
exact retarded spectral moments: 



£<(k,T) = 2n c (k,l>«(k,T) 



1 d 2 n c (k,T) 
! 2 W 2 : 



(69) 



for n = 0, 1, and 2. Since the time derivatives of the 
momentum distribution function are real valued, the last 
term in Eq. I|69() is nonzero only for n — 2. Thus, the 
lesser GKB spectral moments (with n = 0, 1, and 2) are 
equal to the corresponding retarded spectral moments 
multiplied by 2?i c (k, T) plus a term involving the second 
time derivative of the momentum distribution function 
for the case of n = 2. 

Comparison of Eq. (|69|l and the exact expressions for 
the retarded spectral moments [in Eqs. (|19fl . (|23|l and 
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(|26|l ] with the exact expressions for the lesser moments 
[in Eqs. (|37[) and (|38[) ] allows us to conclude that 

the lesser spectral moments for the GKB approximated 
Green's function correspond to an approximation of the 
exact spectral moments by evaluating the operator av- 
erages with a mean-field approximation (the second mo- 
ment contains an additional term with the second time 
derivative of the momentum distribution function). This 
indicates that the GKB approximation will fail as the 
correlations increase, because it produces the wrong mo- 
ments to the spectral functions. Note that this does not 
say the GKB approach is a mean field theory approach, 
it is not. 



IV. CONCLUSIONS 

In this paper, we derived a sequence of spectral mo- 
ment sum rules for the retarded and lesser Green's func- 
tions of the Falicov-Kimball and Hubbard models. Our 
analysis holds in equilibrium and in the noncquilibrium 
case of a spatially uniform, but time-dependent external 
electric field being applied to the system. Our results 
are interesting, because they show there is no average 
time dependence nor electric field dependence to the first 
three moments of the local retarded Green's functions. 
This implies that the value, slope and curvature of the 
local retarded Green's functions, as functions of the rel- 
ative time, do not depend on the average time or the 
field. Such a result extends the well-known result that 
the total spectral weight (zeroth moment) of the Green's 
function is independent of the field. It also implies that 
one will only see deviations of Green's functions from the 
equilibrium results when the relative time becomes large. 
Such an observation is quite useful for quantifying the ac- 
curacy of nonequilibrium calculations. We showed some 
numerical results illustrating this effect for nonequilib- 
rium DMFT calculations in the Falicov-Kimball model. 
The case for the lesser Green's function is more compli- 



cated, and there is both average time dependence and 
field dependence apparent in those moments, although it 
appears that the second moment at half filling may not 
depend on average time. 

We also examined a common approximation employed 
in nonequilibrium calculations, the so-called generalized 
Kadanoff-Baym ansatz. We find the moments in that 
case are similar in form to the exact moments, except 
they involve a mean-field-theory decoupling of correla- 
tion function expectation values when one evaluates the 
operator averages that yield the sum-rule values. This 
implies that the GKB approximation must fail as the 
correlations increase, and such a mean-field decoupling 
becomes inaccurate, because it will have the wrong spec- 
tral moments. 

We hope that use of these spectral moments will be- 
come common in nonequilibrium calculations in order to 
quantify the errors of the calculations for small times. 
We believe they can be quite valuable in checking the fi- 
delity of numerical calculations and of different kinds of 
approximate solutions. 
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APPENDIX A: DERIVATION OF THE SECOND 
SPECTRAL MOMENT FOR THE RETARDED 
GREEN'S FUNCTION 

We show the explicit steps needed to derive Eq. (|25() by 
using Eqs. i|14|) and Q) for the second spectral moment 
of the retarded Green's function. In this case, the second 
spectral moment is equal to 



M 2 R (k,T)= - 



This expression is equivalent to 



/^(k,T) 



dtulm 
dtulm 
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dtulm 

d_ 

dt 



dte 1 ^ 



dtulm 



.dS{t) 



dte lu h 



dt 2 



9{t) ( c k T 



dt 



dte lu}t 5(t) 
dte lU}t i9{t) 



d 



t 



if t c*\ T-- ) ,4( T- - 



t 



T+- ,c 



dt 2 



Mt- 



2<^-c k (T+-),z-c^(T-- 



d 



(Al) 



cu[T+-),i-cl [T-- 



(A2) 
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The first two integrals in Eq. (|A2|) are equal to zero, 
that is, 



cLulx 



at 



= 0, (A3) 



where we replaced derivatives with respect to time by 
commutators with the Hamiltonian. The first term has 
no imaginary part, so it vanishes, as do the second two 
terms, since one can easily show the difference of the two 
operators is Hermitian, and hence has a real expectation 
value. This completes the proof of Eq. I|A3[1 . 



and 



dull 



dte lult S(t) 



c k (T+- T-- 



0.(A4) 



To prove Eq. I|A3[I . we note that the delta- function 
derivative satisfies 



dtf(t) 



dd{t - a) 
at 



-/»■ 



(A5) 



Hence, we can transfer the time derivative of the delta 
function into a time derivative of the other factors in the 
integrand. This time derivative has two terms: (i) the 
first term is the derivative of the exponential, which in- 
troduces an additional factor of ioj (the operator average 
then becomes trivial when we set t = 0, since the anti- 
commutator is equal to 1) and (ii) the second term which 
involves derivatives of the creation and annihilation op- 
erators. Performing the integration over t by using the 
delta function then yields 

- I |^^hmH-Im[({[ Ck (T),W(T)],4(r)} 



|{c k (T),[4(T),W(T)]})] 



(A6) 



To prove l|A4f) . we first perform the integration over t. 
The result is equal to 



Pk(T) 



dulm— 
2 

d_ 



iQfC k (T),ci(T) 



^mi([{[c k (T),W(T)],c[(T)} 



{ Ck (T),[ c t(T),W(T)]}]). (A7) 

The operator in the second line of Eq. I|A7J) is the same 
operator we saw above; it is Hermitian so the imaginary 
part of the statistical average vanishes. This proves i|A4|l . 



To complete the derivation of the second spectral mo- 
ment, we note that only the last integral term in Eq. (|A2|) 
can be nonzero. Substituting the operator time deriva- 
tives by their commutators with the Hamiltonian finally 
gives the expression in Eq. (|25l) of the second spectral 
moment for the retarded Green's function. 
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